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Abstract 

We present a numerical study of the spectrum of the Laplacian in an 
unbounded strip with PT-symmetric boundary conditions. We focus on 
non-Hermitian features of the model reflected in an unusual dependence 
of the eigenvalues below the continuous spectrum on various boundary- 
coupling parameters. 



1 Introduction 

In the last years the theory of quasi-Hermitian, pseudo-Hermitian and VT- 
symmetric operators has developed rapidly, and has been shown to provide a 
huge class of non-Hermitian Hamiltonians with real spectra (c/ the pioneering 
works [24l |4l [23| and the review [3] with many references). Because of these 
recent observations, the condition of self-adjointness of operators representing 
observables in quantum mechanics may seem to be rather an annoying techni- 
cality. However, unless it is met one cannot apply the very powerful machinery 
of spectral theory based on the spectral theorem. 

In particular, one has to restrict to exactly solvable models [22l[27l[25lfT8l[2T] 
or to rely on perturbation and numerical methods to analyse the spectrum 
of the PT-symmetric Hamiltonians. Except for some simple one-dimensional 
examples [2, the majority of the PT-symmctric models studied in the literature 
have purely discrete spectrum. Perturbation methods arc then notably effective 
in determining the dependence of the eigenvalues on various parameters of the 
given model. Although the perturbation approach can even prove that the total 
spectrum is real in some cases [501 ISl IZl IH] : it is limited in its nature and one 
usually has to employ numerical techniques in order to obtain a more complete 
picture of the spectral properties. 

In a way motivated by the lack of a well-developed spectral theory for non- 
self-adjoint operators with non-compact resolvent, in a recent paper [6] Borisov 
and one of the present authors introduced a new two-dimensional T^T-symmetric 
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Hamiltonian with a real continuous spectrum. One of the main questions aris- 
ing within this model is whether the Hamiltonian possesses point spectrum, 
too. Using some singular perturbation techniques adopted from the theory of 
quantum waveguides [T31 H] , the question was given both positive and negative 
answer in [5] , depending on the nature of the effective T^T-symmetric interac- 
tion in a weakly-coupled regime. Moreover, in the case when the point spectrum 
exists, the weakly-coupled eigenvalues emerging from the continuous spectrum 
were shown to be real. We refer to the next Section [5] for a precise statement of 
the spectral results established in [5]. 

In the present paper, we further analyse the question of the existence of 
point spectrum for the Hamiltonian introduced in [5] by numerical methods. 
This enables us to explore quantitative properties of the eigenvalues without the 
restriction to the weakly-coupled regime. The main emphasis is put on peculiar 
characteristics of the model which are related to the non-self-adjointness of the 
underlying Hamiltonian, namely: 

1. highly non-monotone dependence of the eigenvalues on a coupling param- 
eter; as the parameter increases, the eigenvalues emerge from the contin- 
uous spectrum, reach a minimum, sometimes disappear in the continuous 
spectrum, emerge later on again, etc; 

2. broken VT -symmetry; as the coupling parameter increases, the eigenvalues 
may emerge from the continuous spectrum as complex-conjugate pairs, 
collide and become real, move on the real axis, collide again and become 
complex, etc. 

This paper is organized as follows. We begin by recalling the Hamiltonian 
from [5] and summarize the main spectral properties established there. In partic- 
ular, we point out some questions the study of [B] has left open. In Section[3]we 
describe the numerical methods we use. The numerical data arc then presented 
and discussed in Section |4l In the final Section [5] we make some conjectures 
based on the present study. 

2 The model, known results and open questions 

Given a positive number d, we introduce an infinite strip f2 := Mx (0, d). We split 
the variables consistently by writing x = {xi,X2) with xi G M and X2 € (0, c?). 
Let a : R ^ R be a bounded function; occasionally we shall denote by the same 
symbol the function x i— *■ a{xi) on f2. The Hamiltonian Ha we consider in this 
paper acts simply as the Laplacian in the Hilbert space L'^{n), i.e. 

iJ„*:=-A* in n, (1) 

and a non-trivial interaction is introduced by choosing as its domain the set of 
functions 5* from VF^'^(r2) satisfying the following Robin- type boundary condi- 
tions: 

a2*4-iQ!* = on an. (2) 
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Here W'^''^{fl) denotes the Sobolev space consisting of functions on fi which, 
together with all their first and second distributional derivatives, are square 
integrable. As usual, the action of Ha should be understood in the distributional 
sense and © should be understood in the sense of traces [T]. 
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Figure 1: A schematic view of an infinite planar waveguide of width d. The 
Robin conditions ([2]) arc imposed at the boundary. 

Under the additional hypothesis that a possesses a bounded distributional 
derivative, i.e. a G W^'^{M.), it was shown in [B] that Ha is an m-sectorial 
operator satisfying 

Ha = H^a , (3) 

where H* denotes the adjoint of Ha. (If a is merely bounded, it is still possible 
to give a meaning to Ha by using the quadratic-form approach.) Of course. 
Ha is not self-adjoint unless a vanishes identically. However, the relation ^ 
reflects the T^T-symmetry - or, more generally and more precisely, the T-self- 
adjointness - of Ha, with V and T being defined by {V'i>){x) := "^{d — x) and 
T^* := respectively. 

An important property of the operator Ha being m-sectorial is that it is 
closed. Then, in particular, the spectrum cr{Ha) is well defined as the set of 
complex points z such that Ha — z is not bijcctivc. The point spectrum ap{Ha) 
equals the set of points z such that Ha — z is not injective. The continuous 
spectrum adHa) equals the set of points z such that Ha — z is not surjective 
but the range of Ha — z is dense in L'^^fl). Finally, the residual spectrum <T^(Ha) 
equals the set of points z such that Ha — z is injective but the range of Ha — z 
is not dense in L'^{Vl). 

In the following theorem we collect general results about the spectrum of Ha 
established in [B]: 

Theorem 1. Let a e W^^°°{R) and ao G M. Then 

(i) a{Ha) C {z G C : | arg(z)| < 0} with some 6 G [0,7r/2); 

(ii) a, (Ha) = 0; 

(iii) a{Hao) = crdHao) = [^0:°°) where ^iq min{|ao|, Tr/d}; 

(iv) ifa-ao G Co(M), then adHa) = [fil,oo); 

(v) if a Cz Co(K) is an odd function, then <7p{Ha) C M; 
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Here Co(M) denotes the space of continuous functions on K with compact sup- 
port. Note that a G Co(R) n Vl^^'°°(E) imphes that a is Lipschitz continuous; 
conversely, the space of Lipschitz continuous functions on M is embedded in 

The first two properties of Theorem [1] are quite generah (i) holds since Ha 
is sectorial and (ii) is a consequence of the T-sclf-adjointness (we refer to [5] for 
more details). Since the spectral problem for the constant case of (iii) can be 
solved by some sort of "separation of variables" (c/ [51 Sec. 4]), we shall refer 
to it as the unperturbed case; it follows from Theorem [1] that the correspond- 
ing spectrum is purely continuous and positive. As a consequence of (ii), (iv) 
and (v) , we get a result about the reality of the total spectrum in the perturbed 
case: 

Corollary 1. Let a £ Co(R) n M^i'°°(R) he an odd function. Then 



The result stated in part (iv) of Thcorem[T]makes rigorous the heuristic state- 
ment that "the continuous spectrum depends on the properties of a Hamiltonian 
interaction at infinity only" . On the other hand, it is well known - and this 
already for one-dimensional self-adjoint models [17] - that the point spectrum 
may be highly unstable under a perturbation of an operator with non-compact 
resolvent. In [6], the point spectrum of Ha was analysed perturbatively in the 
weakly-coupled regime: 



where /3 is assumed to be real-valued and e plays the role of the small param- 
eter. Here Cq{M.) denotes the space of functions on K which, together with all 
their first and second derivatives, are continuous and have compact support. 
The main interest was focused on the existence and asymptotic behavior of the 
eigenvalues emerging from the threshold /ig of the continuous spectrum due to 
the perturbation of Hao by ef3. 

Before stating the main results of [5] about the weakly-coupled eigenval- 
ues, we need to introduce some notation. Recalling the definition of /ig from 
Theorem [1] (iii), we next define 



a{Ha) C M. 



a = uq + e P , 



with ao G M , e > , /? G C^(M) , 



(4) 



^1 := max{|Q!o|, Tr/d} and fij :~ irj/d for j > 2. 
To these numbers we associate a family of functions {■i/'jlf^o 

ipj{x2) ■■= cos{fijX2) - i — sin(^jX2) . 
Let {vj}°^Q be a sequence of auxiliary functions given by 



(5) 
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Finally, denoting (/) = /(xi) dxi for any / G L-'^(R), we introduce a con- 
stant T, depending on /3, d and ag, by 



2ao7r2coti^§^,^ ^ 8^2 ^ fih{(3v2j) 



2ao(/3wo) + -1- 2^ 4 2 tan if |ao| < -r 



(6) 



Now we are in a position to take over from [S] : 
Theorem 2. Lei a &e ^iuen by 

(A) //ao = 0, then Ha has no eigenvalues converging to /Iq as e ^ 0. 

(B) Let < |ao| < 7r/d. 

1. If aQ{P) < 0, i/ien i/iere exists the unique eigenvalue Xg of Ha con- 
verging to as e 0. This eigenvalue is simple and real, and 
satisfies the asymptotic formula 

Xe^fil^ e^alip)^ + 2e^aoT{p) + 0(6^) . (7) 

2. If ao{j3) > 0, then Ha has no eigenvalues converging to /j,q as £ — > 0. 

3. // = and r > 0, then there exists the unique eigenvalue of 
Ha converging to as e — > 0. This eigenvalue is simple and real, 
and satisfies the asymptotics 

Xe= ^il-e\^' +0{e^). (8) 

4. // (/3) ~ and r < 0, then Ha has no eigenvalues converging to Hq 
as e — > 0. 

(C) Let \ao\ > n/d and aQd/i: ^ Z. 

1. If T > 0, then there exists the unique eigenvalue X^ of Ha converging 
to /ig as e ^ 0, it is simple and real, and satisfies the asymptotics ([5]). 

2. If T < 0, then Ha has no eigenvalues converging to /ig as £ ^ 0. 

The method of [5] gives also the asymptotic expansion of the eigenfunctions 
corresponding to the weakly-coupled eigenvalues: 

Theorem 3. The eigenfunction corresponding to any eigenvalue X^ from 
Theorem\^ can be chosen so that it satisfies the asymptotics 

^,{x)^Mx2) + 0{e) (9) 



W^'^{fl n {a; : < a}) for each a > 0, and behaves at infinity as 



*e(a;) = exp-V/^o-^-l^il ^o(a;2) 4-C'(exp-VA'o--^.l^il) , ^ +00 . (10) 
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Theorems [TH2] summarizing the spectral analysis performed in [5] leave open 
the following particular questions: 

(Ql) Can the cases (B4) and (C2) of Theorem [2] occur? That is, can the con- 
stant r be negative for a certain combination of d, ao and (31 (Sufficient 
conditions for the positivity of t exist [SI Props. 2.1-2.2].) 

(Q2) What happens in the case (C) of Theorem [5] if the condition aQdjn ^ Z 
is not satisfied? Is it just a technical hypothesis? 

(Q3) Is there any point spectrum in the case of Corollary [T]/ 

(Q4) What is the dependence of the weakly-coupled eigenvalues of Theorem [5] 
as the parameter e increases? 

(Q5) Do the eigenvalues remain real for large el 

(Q6) Can one have more eigenvalues? Can they be degenerate? What is the 
dependence of the number of eigenvalues on e? 

(Q7) Are there any eigenvalues emerging from other thresholds j > 1? Can 
they emerge from other points of the continuous spectrum, different from 
the thresholds fJ,^, j > 0? 

The main goal of the present paper is to provide answers to some of these 
questions by a numerical study of the spectral problem. 

3 Numerical methods 

In order to get the dependence of the bound states on parameters like e,d,ao, 
etc, numerically we used two independent methods. When a is a simple step- 
like function {e.g. symmetric or asymmetric square well), we treat the problem 
by mode matching method. It takes into account the asymptotic behaviour of 
solution explicitly and can serve thus as a useful check when we apply the other 
method, viz, the spectral method. This method is more robust and we use it 
for more general a. We arrived at an excellent agreement in cases when both 
methods are applicable. 

3.1 Mode matching method 

Let us begin with mode matching. The most general situation we want to 
describe is shown in Figure O Fix negative and positive numbers L_ and L+, 
respectively. In the asymptotic regions, i.e. xi < L_ and < xi, we assume 
a{xi) = ao, while in the central parts we have a{xi) = a_ if L_ < xi < and 
a{xi) = a+ if < xi < L+. 

Let {nf}""^^ denote the sequence of numbers {nj}j°,Q with ao being replaced 
by a±. In the same way we define the sequence of functions {ipfyj^Q by replac- 
ing ao by a± in ([5]). In order to make the notation more consistent, hereafter 
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Figure 2: The mode matching approach. A particular Ansatz ([TT|) for an eigen- 
function ^! of Ha corresponding to A is chosen in each subregion and the smooth 
matching (fT2|l is required at the boundaries separating the subregions. 



we write /i^ and t/j^ instead of fij and i/jj , respectively, and introduce a common 
index l S {0, +, — }. 

In each of the regions where a is constant, the spectral problem — A^" ~ A^*, 
with ^ satisfying the required boundary conditions, can be solved explicitly [SI 
Sec. 4] by expanding ^ into the "transverse basis" {V'^lj^c where l depends on 
the region. More specifically, we use the following Ansatz for an eigenfunction ^ 
of Ha corresponding to A: 



dj eVK')'-^-i ^0(x2) if XI G (-00, L_) , 



if xi e {L_,Q) , 
if xie(0,L+), 



ififixi) ■.= coa{JX-{^if)^Xl)+B± sin(JA- 



OO 



(11) 



^a,e-V("")'-^-i^0(a;2) if xie(L+,+oo) 
i=o 



where 



Standard elliptic regularity theory implies that any weak solution '5 to 
— A^P = A^" is necessarily infinitely smooth in the interior of fl. In particu- 
lar, we must match the functions from the Ansatz smoothly at xi = i_, 0, L+, 
i.e. we require 



vI'0(L_,a;2) = *-(L_,X2), 
*-(0,a:2) = *+(0,X2), 

Vl/+(L+,.T2)=*°(i+,X2), 



ai*?(i_,2;2) -9i*,-(L-,X2), 

ai*-(0,2;2) -9i*+(0,X2), (12) 
ai*+(L+,X2) =ai*5t(L+,a;2), 



for every X2 G (0, d). 

If {"(Ajlj^o formed an orthonormal family, the next step would consist in 
employing the orthonormality and reducing (|12p into a system of algebraic 
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equations for the coefficients aj,bj,Cj,dj. However, since the family {^pj}'^^ 
is actually formed by eigenfunctions of a transverse eigenvalue problem which is 
not Hermitian (unless =0), it is clear that the functions tpj are not mutually 
orthogonal in general. Instead, we use the property that {4'j}j^o and {4'j}j^o 
form a complete biorthonormal pair 
eigenfunctions of the adjoint transverse problem: 



, where (j)^ are properly normalized 



The normalization constants can be chosen as follows 
2i«. 2(/ii)2 



^" ■ 1 - exp (-2ia,d) 



where j ^ 2, (jo,ji) = (0,1) if |aj < 7r/d and (jo,ji) = (1,0) if > n/d 
(if at = 0, the fraction in the definition of A^^^ should be understood as the 
expression obtained after taking the limit ^ 0). Then, in particular, we have 

Vz,jeN, (0,J,V';O = ^y, (13) 

where (•, •) denotes the inner product in L^((0, c?)), antilinear in the first factor 
and linear in the second one. 

Now, multiplying by 0°, integrating over a;2 € (0, d) and employing 
in the asymptotic regions, we can eliminate the coefficients Uj and dj by means 
of the relations 



3=0 



6,^+(i+)(</>^,^+), 



(14) 



^c,^-(i_)(^°,V'-), 



for every i e N, and reduce thus the number of conditions to be fulfilled. We 
finally arrive at an infinite-dimensional homogeneous system 



/ mil 

"1-21 



V 



mi2 



77142 





17123 
77144 



\ 


( ' 


\ 






\ 




c 













bB+ 











J 


\ cB^ 


) 




I 





(15) 



Here 6, c denote the infinite vectors formed by bj , Cj , respectively, and the sub- 
matrices 771^1, are given by 



mil 
mi2 



m33 := y^A-(A.+ )2(0°>;), 



'T134 
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TO21 



I^^A^W^ COS (L+y^A-(/i+)2) 

- ^X-i^y sin (L+y^A-(^+)2)^ (0^, ^+ 



77142 



77144 



m23 - A sin (L+^A - 

+ ^A - cos (L+y'A-(M+)2)) (0°, V^+) , 

- ^ cos (L_ y^A-(^T)2) 
+ ^A-(m7)2 sin (i_ ^A^O^)) (0°, V") , 
y^(A.?)2-Asin(L_y^A-(A.7)2) 

- v/A-(/i7)2 cos (i_y'A-(M7)2)) (^°, , 

where the right hand sides should be understood as the infinite matrices formed 
by the respective coefficients for i,j g N. Our numerical approximation then 
consists in approximating the infinite system by using finite submatrices for 
i,j € {0, . . . , A^} with N large enough. 

In order to have a nontrivial solution we require 

/ 77111 "^12 \ 

det "1^1 7.23 ^ (ig) 

77133 ^34 ^ ' 

\ 77142 70,44 / 

which gives an implicit equation for A as the unknown. Having found A, we can 
then calculate the coefficients aj,bj,Cj,dj,B±j from and (|14p . 

If a+ = a_ and L+ = — i.e. a is a symmetric square well, then ([T5| can 
be reduced to 

77121 + 77142 



71123 - 'TI44 

The solutions have different symmetry with respect to xi —xi, the even 
solutions are formed only of cosines, the odd of sines. 

Remark. In principle, it is possible to extend the present method to an ar- 
bitrary piece-wise constant function a, provided that the number of matching 
interfaces is finite. On the other hand, the more matching conditions the bigger 
size of the matrix of (jlSp and thus the higher (numerical) price one must pay. 

3.2 Spectral method 

In order to treat the waveguide with a general a in the boundary conditions, 
we decided to use spectral collocation methods. They provide a reliable and 
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rapidly converging tool easily applicable to our model. A very useful software 
suite has already been published [26] . It could be adapted to this problem. We 
approximate the operator by a series of operators defined on a finite domain 
[a;™", a;™""^] x [0,d] with Dirichlct boundary conditions on {x"""} x [0,d] and 



First, we can form the differentiation matrices in each variable separately and 
then combine them to the two-dimensional problem. The infinite domain in xi- 
variable suggests that we could approximate it by grid points chosen as the roots 
of Hermite polynomials and as an interpolant we take Lagrange polynomial. 
There is an additional parameter (the real line can be mapped to itself by a 
change of variable xi = > 0), which can be used to optimize the choice 

of the grid points together with variation of the number of roots Ni . The roots 
span the interval [^i,^ArJ, — ^^i = £.Nii cluster around the origin, and grow as 



Another possibility is to use Fourier differencing. Wc form a uniform grid in 
^_^7nax ^ ^maxj g^^^ sincc the solutions decay exponentially, we can theoretically 
extend it periodically across this interval to the whole R. The interpolant is a 
trigonometric function. 

In both approaches the interpolant is an infinitely differentiable function. 
Deriving it and taking the derivatives in the grid points we get the differentiation 
matrices. Imposing the homogeneous Dirichlet boundary conditions consists in 
deleting the first and the last rows and columns of the differentiation matrices. 

The transversal variable is confined to a finite interval [0, d] and it is possible 
to scale it to [—1, 1]. To implement the boundary conditions we prefer to incor- 
porate them into the interpolant. It requires to use Hermite interpolation, which 
takes into account derivative values in addition to function values. The use of 
the roots of Chebyshev polynomials rjk = cos((fc — 1)tt/ {N2 — 1)), k = 1, . . . , N2 
as the grid points is common here. For details wc refer the reader to p6] . 

Now, it remains to form differentiation matrices that correspond to partial 
derivatives entering the Laplacian. Having set up a grid in each direction we 
combine them into the tensor product grid. Then a closer inspection shows that 
df D*^^' where I is an N2 x N2 identity matrix and 9| is constructed in 

a similar way (it is necessary to take into account that the boundary conditions 
change with xi). 

Applying the spectral discretization we converted the search for eigenvalues 
of Ha to a matrix eigenvalue problem. 

4 Discussion of numerical results 

Existence of eigenvalues below the threshold of the continuous spectrum and 
their behaviour for weak perturbations was already proved in [5] . Our calcula- 
tions confirm it and demonstrate that the spectrum of eigenvalues is consider- 
ably richer. 

A typical dependence of an eigenvalue on the perturbation parameter e is 
shown in FigurcO Here we perturbed a waveguide of width d ~ 2 and ao = 1/3 




) for iVi 00. 
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Figure 3: Comparison of dependence of eigenvalues on e for T^T-synimetric 
and self-adjoint waveguides. The left figure shows the T'T-symmetric case with 
a(x\) = 1/3— e exp(— x^). Here the blue (respectively magenta) curve represents 
the eigenvalue (respectively the asymptotic formula ([7]) up to the e'^-term). The 
right figure shows the eigencurves in the corresponding self-adjoint situation, 
obtained by replacing ia i— > a in . c? = 2 in both cases. 

by a Gaussian shape, i.e. we took /3(a;i) = — exp(— x^) in ([4]). Since < ag < 
Tr/d and (/3) < 0, we deal with the case (Bl) of Theorem [2l We observe that 
the asymptotic expansion ([7]) is fairly good. It is striking, however, that the 
dependence of the eigenvalue on e is highly non-monotonic: The eigenvalue 
appears at some value of e (in this case it is e = 0), reaches a minimum, and 
then returns to the continuous spectrum. We found such a behaviour in all 
cases we studied, viz^ various shapes of symmetric and asymmetric wells, and 
Gaussians times polynomials. This provides an interesting answer to (Q4) from 
the end of Section [51 

It is worth noting that this behaviour differs from that in the self-adjoint 
waveguide obtained simply by omitting the imaginary unit in As shows 
the second graph in Figure [3l in the self-adjoint case all the energy levels are 
increasingly more bound when e increases. 

On the other hand, we checked that the eigenvalues are decreasing as func- 
tions of L := ±L± for the symmetric square-well profile a+ = a_ of Section [3.1l 
in the regime < a± < < n/d. This is reasonable to expect since as L — s- cx) 
the eigenvalues should approach (a±)^, i.e. the threshold of the continuous 
spectrum of the unperturbed waveguide H^^ . 

An answer to questions from (Q6) is provided by FigurclH It corresponds to 
the critical case (B3) of Theorem [2] with P{xi) ~ —{xf + bxi — 5) cxp(— x^/lO), 
ao = V2, and d = 2; the parameter b changes the asymmetry of (3. In addition 
to the weakly-coupled eigenvalue of Theorem [21 there are also other eigenvalues 
emerging from the continuous spectrum as e increases. The lower figure (case 
b 3.25) shows that there might be eigenvalues existing in disjunct intervals 
of £ (the red curve) . By diminishing aQ we can achieve the situation when there 
is only one eigenvalue with a similar behaviour: it emerges from the threshold of 
continuous spectrum (at e = 0), reaches a minimum, returns to the continuum. 
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Figure 4: Dependence of eigenvalues on e in the critical case (/?) = 0, r > 0. 
Here a{xx) = \/2 — e{x\ + hxi — 5) exp(— a-f /lO) and d = 2. The upper figure 
corresponds to 6 = 1.5, the lower one ioh — 3.25. All the crossings are avoided. 



reappears later on, and returns finally to the continuum. 

Another typical feature is that the energy levels do not cross. We saw al- 
ways avoided crossings (at least in the unbroken ■pT-regime), i.e. the order of 
levels remains unchanged and the non-monotonicity of the excited eigenvalues 
is preserved. 




Figure 5: Broken T'T-symmetry. The left figure shows the dependence of eigen- 
values on £ in the case a{xi) = 7r/2 — eexp(— x^/10) and d = 2. Here the 
dotted line is used to plot the real part of the eigenvalues if they form complex 
conjugate pairs instead of being real. The right figure shows the trajectory of 
a pair of (complex) eigenvalues in the complex plane. Here the thick black line 
marks the continuous spectrum. The pairs of dots show positions of eigenvalues 
for different values of e: 0.1 (blue). 1.3 (red), 2.1 (green), 2.5 (black), and 8 
(magenta). An animation can be found at the website [19j . 
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The spectrum in Figure is remarkable from two points of view. First, it 
corresponds to the case of (Q2) from the end of Section [2l since the constant 
ao is chosen so that its square coincides with the threshold of the continuous 
spectrum, which is (7r/2)^ « 2.47 for d — 2. Second, we see that this situation 
provides a negative answer to (Q5), i.e. the PT-symmetry can be broken, and a 
partial answer to (Q7). Let us comment on the behaviour depicted by the cyan 
curve. There is a critical value of the parameter e for which there emerges a pair 
of complex conjugate eigenvalues from the continuous spectrum (we suspect that 
they emerge due to a collision of two embedded eigenvalues). As e increases, 
the eigenvalues propagate in the complex plane (this is indicated by the curves 
joining the blue and red dots in the right figure; the dotted curve in the left 
picture traces the common real parts of the eigencurves) till they collide on the 
real axis and become real. Then they move on the real axis (as the green dots) 
in opposite directions till they reach turning points (each of them for different 
value of e), starts to approach each other, coalesce again and continue as a pair of 
complex conjugate eigenvalues (indicated by the black and magenta dots) until 
they disappear in the continuous spectrum. The behaviour of the eigenvalues 
depicted by the red curve in the left picture is more difficult in that one of 
them seems to have the turning point inside the continuous spectrum. Because 
of the collisions we see that the eigenvalues can actually be degenerate if the 
■PT-symmetry is broken, providing a positive answer to one of the questions 
from (Q6). 

Let us mention that the behaviour of the eigenvalues in the regime of broken 
T'T-symmetry exhibits certain similarities with the over-damped phenomena 
as regards the spectrum of the infinitesimal generator of the semigroup associ- 
ated with the damped wave equation [101 HH IE] • This indicates the unifying 
framework of Krein spaces behind these two problems [101 lU • 



5KW 3(V) 




"1 



Figure 6: Real and imaginary parts of an eigenfunctions corresponding to the 
smallest (positive) eigenvalue of Ha for a{xi) ~ 1/3 — 0.65 exp(— 0.025xf ). 

In Figure [6l we present an example of eigenfunction corresponding to the 
case (Bl) of Theorem [2l We check that the behaviour of the eigenfunction is in 
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perfect agreement with the asymptotic results of Theorem [3] 

Even if the corresponding eigenenergies are real, the non-Hermiticity pre- 
vents from choosing the eigcnfunctions real. Since the latter is in particular 
true for the lowest eigenvalue, it does not make sense to speak about the super- 
and sub-harmonic properties of the corresponding eigenfunction (which hold in 
the self-adjoint case). However, although there is no variational characteriza- 
tion of eigenvalues in the present model, numerically wc observe that the real 
and imaginary parts of the eigenfunction corresponding to the lowest eigenvalue 
are super-harmonic separately in the regime of unbroken T^T-symmetry. This 
follows, of course, from the observations that they do not change sign and that 
the spectrum is positive. 
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Figure 7: Dependence of r on parameters defining a in the step-like situation 
of Section 13.11 with d = 2. The first two figures deal with an antisymmetric 
square well in the regime |q!o| < 7r/rf, while the last one deals with a symmetric 
square well in the regime oq > n/d. The first (respectively second) figure shows 
the dependence of r on the width ±L± =: L (respectively on the coupling a) 
for fixed uq ~ 1/3 and a± = ao T 1 (respectively fixed ao — 1, ±i± 2 and 
variable a± = ao T (5). In the last figure we show the dependence of r on ao for 
L = 10 and variable a± = (ao — 1); here the dotted line corresponds to ir/d. 

Finally, in Figure[7]we visualize the dependence of the complicated quantity r 
defined in ([6]) on various parameters. In particular, we see that it changes sign, 
giving a positive answer to (Ql). Consequently, all the cases of Theorem [5] for 
the critical case = and for the regime |ao| > T^/d can be achieved. We 
also see that the first figure is in qualitative agreement with an analytic result 
of P Prop. 2.1]. 



5 Conclusion 

In this paper we tried to enlarge our knowledge of the point spectrum of a non- 
Hermitian T'T-symmetric operator introduced in [6] by analyzing it numerically. 
We confirmed theoretical results obtained in [6j by perturbation methods, and 
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showed that they actually hold under much milder conditions about a. 

Besides this, it turned out that the operator can model a fairly wide range of 
situations. Indeed, its spectrum is very rich, and certain properties we found are 
unusual when we compare them with the standard self-adjoint cases. Among 
them we would like to point out the non-monotonic dependence on the strength 
of perturbation and the existence of the regime of broken PT-symmetry. We 
hope that this study will stimulate further theoretical endeavour to extract and 
prove the salient features. 

In particular, based on the present numerical analysis, we conjecture that 
there will be no other spectrum except for the continuous one if the parameter e 
is sufficiently large. At the same time, we were not able to find any discrete 
eigenvalues in the case of Corollary [U i.e. (Q3) from the end of Section [2] seems 
to have a negative answer; the statement (v) of Theorem [1] would be trivial, 
then. Our numerical experiments also indicate that the condition mentioning 
in (Q2) is indeed just a technical hypothesis in Theorem [2lC, in the sense that 
it does not influence the existence/non-existence of weakly-coupled eigenvalues. 

More generally, the existence of eigenvalues in the present model seems to 
have a nice heuristic explanation. We observe that the discrete spectrum be- 
haves in many respects as that of a one-dimensional Schrodinger operator gov- 
erned by the first-transverse-eigenvalue potential, i.e., — A -I- min{a^, Tr^/d^} in 
L^(R). Of course, this self-adjoint idealization is just approximative and can- 
not explain, in particular, the existence of non-real eigenvalues. However, it 
provides an insight into the non-monotonicity behaviour, the absence of point 
spectrum for large e, the positivity of (the real part of) the spectrum, etc. It also 
formally explains the condition from the statement 1 (respectively 2) of Theo- 
rem[2jB, since this actually implies that the potential is attractive (respectively 
repulsive). 

In this paper we were mainly interested in the eigenvalues emerging from 
the threshold of the continuous spectrum. A complete answer to the first 
question of (Q7) can be provided by a perturbation method similar to that 
of [B]. However, a more detailed analysis of the continuous spectrum would 
be still desirable. In particular. Figure S] suggests that there can be embedded 
eigenvalues for larger values of the coupling parameter. 

Finally, let us point out that the question of a direct physical motivation 
for the Hamiltonian Ha remains open. In this paper we were rather interested 
in consequences of the non-self-adjointness on spectral properties of this spe- 
cific model in the context of T^T-symmetric quantum mechanics. On the other 
hand, motivated by problems in semiconductor physics, similar self-adjoint, re- 
spectively non-self-adjoint but dissipative, Robin-type boundary conditions has 
been considered recently in [T5], respectively in [TH]. In a different context, the 
present (T'T-symmetric) Robin-type boundary conditions imply that we actu- 
ally deal with the Helmholtz equation in an electromagnetic waveguide with 
radiation/dissipative boundary conditions. 
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